SOME GRADIENT ESTIMATES FOR THE HEAT EQUATION ON 
DOMAINS AND FOR AN EQUATION BY PERELMAN 



QI S. ZHANG 

Abstract. In the first part, we derive a sharp gradient estimate for the log of Dirichlet 
heat kernel and Poisson heat kernel on domains, and a sharpened local Li-Yau gradient 
estimate. 

In the second part, without explicit curvature assumptions, we prove a global upper 
bound for the fundamental solution of an equation introduced by G. Perelman, i.e. the 
heat equation of the conformal Laplacian under backward Ricci flow. Further, under 
nonnegative Ricci curvature assumption, we prove a qualitatively sharp, global Gaussian 
upper bound. The idea is to combine the Nash and Davies heat kernel estimate with a 
Sobolev imbedding by Hebey, together with a Hamilton type gradient estimate. 
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1. Introduction 

The goal of the paper is to establish certain new point-wise or gradient estimates for 
the heat equation in both the fixed metric and the Ricci flow case. Gradient estimates for 
the heat equation are important tools in geometric analysis as pointed out in the papers 
[LY], [H], [CH], [ACDH] and others. In this paper, for the fixed metric case, two gradient 
estimates are proven. One is a sharp gradient estimate for the log of Dirichlet heat kernel 
and Poisson heat kernel on domains. This can be viewed as a boundary version of the well 
known Li-Yau gradient estimate (see Section 2 for a restatement). It can also be viewed 
as another step in the long running process of heat kernel or Poisson kernel estimate 
starting with the Gaussian formula and Poisson formula. As far as boundary gradient 
estimate is concerned, only the Neumann boundary case was treated in [LY] and [Wa], 
The Dirichlet case is different in that solutions vanish on the boundary. Therefore, the 
estimate is different from both the Li-Yau theorem and its generalization in [Wa]. The 
result seems to be new even for Euclidean domains. The other result is a sharpened local 
Li-Yau gradient estimate that matches the global one. As well known, even for manifolds 
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with nonnegative Ricci curvature, the local Li-Yau estimate differs with the global one by 
a nontrivial factor. Here we show that this factor can be chosen as one. Hence the global 
and local Li-Yau estimate are identical. We expect the result to have applications in heat 
kernel estimate on manifolds. These two results are presented in section 2. 

In sections 3 and later we turn to the case when the metric evolves by the Ricci flow. 
In the fundamental paper [P], Perelman discovered a monotonicity formula for equation 
(4.0) below, which can be regarded as the heat equation of the conformal Laplacian under 
backward Ricci flow. Perelman's formula can be thought of as a gradient estimate. Using 
this estimate together with the reduced distance function, he then applies the maximum 
principle to prove a lower bound for the fundamental solution of (4.0), whenever it exists. 
The reduced distance incorporates the scalar curvature as an integral part. However it 
does not seem that the maximum principle alone will yield a two sided bound. Here we 
will use the Nash method to show that certain long time upper bound also holds. The 
bounds involve more classical geometric quantities such as the best constants in Sobolev 
imbedding, which depend only on Ricci curvature lower bound and injectivity radius. 
Under no explicit curvature assumption, we prove a global on-diagonal upper bound for 
the fundamental solution on (4.0). The bound is good in the sense that it matches the on 
diagonal bound in the fixed metric case. However, we are not able to prove a good off- 
diagonal bound without further assumptions on curvature. Nevertheless, the on-diagonal 
upper bound does not have the usual, trouble making exponentially growing term even 
when the Ricci curvature changes sign. When the Ricci curvature is nonnegative, we 
obtain a qualitatively sharp Gaussian upper bound. This is presented in Section 5. 

We will use the following notations throughout the paper. M denotes a compact Rie- 
mannian manifold without boundary, unless stated otherwise; g, Rij will be the metric 
and Ricci curvature; V, A the corresponding gradient and Laplace-Beltrami operator; c 
with or without index denote generic positive constant that may change from line to line. 
In case the metric g(t) evolves with time, then d(x, y, t) will denote the corresponding 
distance function; d/j,(x,t) denotes the volume element under g(t); We will still use V, A 
the corresponding gradient and Laplace-Beltrami operator, when no confusion arises. 

To close the introduction, we point out that all results in the paper are stated for com- 
pact manifolds or bounded domains. However similar results can be proven for noncompact 
manifolds under appropriate assumptions near infinity. 



2. LOG DERIVATIVE ESTIMATES FOR DlRICHLET HEAT KERNEL AND POISSON HEAT 

KERNEL ON DOMAINS. 

There have been several log gradient estimates available for the heat kernel on complete 
manifolds, compact manifolds without boundary and for the Neumann heat kernel. In the 
introduction, we mentioned the papers [LY] and [Wa]. For compact manifolds without 
boundary, we refer the reader to the papers [Sh], [H] (Corollary 1.3), [No], [MS], [Hs] and 
[ST]. However, an estimate for the Dirichlet heat kernel is clearly missing. This is done in 
the next theorem. The estimate is sharp in general as can be seen from the heat kernel 
formula for the Euclidean half space. Let us mention that for bounded domains, the large 
time behavior of heat kernels is determined by the first eigenvalue and eigenfunction. So 
we will only deal with the most interesting, small time case. 
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Theorem 2.1. Let D be a bounded C 2 domain in a Riemainnian manifold and G = 
G(x, t; y, 0) and P = P(x, t; y, 0) be the Dirichlet heat kernel and Poisson heat kernel 
respectively. Also let p(x) = dist(x, dD) and d(x, y) be the Riemannian distance. Given 
T > 0, there exists a constant C depending on T and D such that 



(2.1) \V x logG(x,t;y,0)\< 



+ P (x)>Vt; 



p(x) • 

Vt ^ 1 ~7T 

for all x,y £ D and < t < T; and for all x £ D, y £ dD and < t < T, 
(2.2) |V.logP(x, t ; 9 ,0,|<{# | ' i + ,g ) ] £ ^ w>v/J 



Proof. 



Let us prove (2.1) first. The proof of (2.2) is similar and will be sketched later. As 
explained earlier, the most interesting case for the derivative estimate is for small time. 
Hence we can take T to be sufficiently small. Here we will take T so small that the 
boundary Harnack principle of [FGS] holds when p(x) < 2T. Here we notice that even 
though the boundary Harnack principle was proven in the Euclidean case in that paper, it 
is still valid in the current case. This is so because we can cover the boundary of D by a 
finite number of metric balls with radius less than the injectivity radius. And then we can 
convert the Laplace-Beltrami operator into an elliptic equation with smooth coefficients 
in R n . 

For a fixed to £ (0, T) and y £ D, we write 

f(x,t) = G(x,t;y,-to), x£D,t>0; 

(2.3) n t0 = {(z, t) I x £ D, < r < t Q , p{z) > ^ }. 

Fixing (x, t) £ D x [0, to] — ^t -> we can a PPly the gradient estimate in Theorem 1.1 of [SZ] 
on the cube 

Q x , t = B{x,p{x)) x[t- ?4^-,t] cDx [-t , to]. 

This gives us 

(24) ^<JL (l + log ^U 

f(x,t) p(x) y f{x,ty 

Here A = supQ xt f. For a proof of (2.4) and that of Theorem 1.1 in [SZ], please go to 
Theorem 3.1 in the next section, which contains Theorem 1.1 in [SZ] as a special case. 
Now we apply the standard Harnack inequality of [LY] on manifold to reach 

A = sup Qxt f < cif(x, t + p{x) 2 ). 

Then the boundary Harnack inequality of [FGS] gives us 

f{x,t + p{x) 2 ) <c 2 f(x,t) 

since / vanishes on 3D x (—to, to)- Therefore 

(2.5) A = su PQxt f <c 3 f(x,t) 
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Substituting (2.5) to (2.4), we deduce, for (x,t) G D x [0,t ] - Q to , 
(2.6) 



\Vf(x,t)\ < C 



/(rr,t) " p(x) 

This proves the first part of (2.1). 

Next we work in Qt - Let us observe that on the sides of dQt , there holds p(x) = \[t. 
Hence for such x, inequality (2.6) becomes 

|V/(M)| < C 



i.e. 



(2.7) 



f( X ,t) - y/t 



\Vf(x,t)\ 2 C r- 



Let m = sup n /, then for any b > 0, we have 

/ log ^j- > / log b. 

Now we use the calculation in the proof of Theorem 1.1 in [H] (pll5) to reach 

A(/log^)-Wlog^) 



(2.9) 



(Af -d t f)logb + (A- d t )(f log j 
|V/| 2 



Also 



(2.10) 



|V/| 2 



difdjf 



f 



+ 2Ri 



djfdjf 
f 



> -2K 



f 



Here — K is the lower bound of the Ricci curvature. Therefore, for 

t |V/| 2 6m 

h= i+2Kt f ~ /log r 

we have 

A/i(x, t) - 9 4 /i(x, i) > 0, (x, t) G fi xo . 

We t = 0, it is clear that h < 0. On the sides of <9f^ , i.e., when p(x) = \/t, one can choose 
b sufficiently large so that 

Here we just used (2.7) and the constant C is from there too. Therefore h < on the 
sides of dQ.t . When T is sufficiently small, we know that dflt is connected and we can 
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apply the maximum principle on this time dependent domain to conclude that h(x) < 
in Q to . i.e. 



f2.ll) 



t Ef <c ( r,*)/i„ g ^ 



In particular, this shows, with 

t = t ,f = G(x, t ; y, -t ) = G(x, 2t ;y, 0), 

that 



V x G(x,2t ;y,0)\ < C(T,K) 



G(x,2t ;y,0) 
when p{x) > \f2t$. Here 



l + 4/lo; 



m 



G(x,2t ;y,0) 



m = sup / 



sup 



f(z,T) 



Sko p(z)> v / 7,0<r<to p{z)>^Ffl<T<t a 



sup G(z,T;y,-t ). 



Making a change of variables to — * to/2, we have 



(2.12) 



\V x G(x,t ;y,0 
G(x,t ;y,0) 



< 



C{T,K) 



I + 4/I0; 



m 



G(x,t ;y,0) 



when p(x) > \/io. 

By the Dirichlet heat kernel upper bound in Davies [Da], we know that 



(2.13) m < C. . , , . , r 

K J ~ \Vh )\B{y,JTo)\ 

By the lower bound estimate in [Z], there holds 



(2.14) 



G(x,t ;y,0) > Ci 



-cd(x,y) 2 /t 



K^M )\B{y,*M)\ 

Here we note that the lower bound was proven under the assumption that the Ricci 
curvature is nonnegative. However for short time behavior this assumption is not necessary. 
Substituting (2.13) and (2.14) to (2.12), we obtain 

d(x,y)~ 



(2.15) 



|V x G(x,*o;y,0)| < C(T,K) 



1 + 



Vto~ . 



G(x,t o ;y,0) y/t 

when p(x) > Now (2.1) follows from (2.7) and (2.15). 

To prove (2.2), let us recall the results in [Da] (upper bound) and [Z] (lower bound): 
there exists c\ and C2 such that 



1 / p(x) 



ci V y/t 



A 1 



\Vt 



A 1 



1 



\B(y,Vt)\ 



-d(x,y) 2 /(c 2 t) 



< 



G(M;s ,„ )£ci (*) A1 )(^ A1 )_i 



Vt)\ 



-c 2 d(x,y) 2 /t 



for all x, y G D and < t < T. 

Given y G dD, the Poisson heat kernel is defined as 

d 

P(x, t; y, s) = -—G(x, t; y, 0). 
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Therefore one has the two-sided bound 

-Al) , , 1 =- e -^,y) 2 /(c2t) < />(,../■„()) ; r ,i: _ A ;i j_ _, 



ciVVt J\B(y,Vi)\ ~ ^ ' ^ ~ V V* J\B(y,Vi)\ 

The rest of the proof for (2.2) is identical to that of (2.1). □ 

Our next theorem provides a sharpened local Li-Yau estimate. In 1986 Li and Yau 
proved the following famous estimate. 

Theorem (Li- Yau [LY]). Let M be a complete manifold with dimension n > 2, Ricci{NL) > 
—K, K > 0. Suppose u is any positive solution to the heat equation in B(xq,R) x [to — 
T, to] C M x [to — T, to]- Then, for any a G (0, 1), there exists a constant c = c(n, a) such 
that 

|Vu| 2 u t c c 
a o < + tt; + cK i 



V? u 



R 2 T 



inB(x ,R/2) x [t -T/2,t ]. 

Moreover, ifM. has nonnegative Ricci curvature and R = oo, i.e. B(xo,R) = M, then 

|V^| 2 _ Ut < Cn 

u 2 u - T' 



Let us observe that, even in the case of nonnegative Ricci curvature, the first local 
estimate does not match the second global estimate completely, due to the presence of the 
parameter a < 1. Here we show that a can be taken as 1 modulo a lower order term. We 
mention that our estimate in the next theorem is new only in the local sense. The global 
estimate was already proven in [Y] by a using a more involved quantity. The very short 
proof, simpler than previous ones, is based on a modification of an idea in [H] and the 
cut-off method in [LY]. 

Theorem 2.2. Let B(xq,R) be a geodesic ball in a Riemannian manifold M with dimen- 
sion n > 2 such that Ricci\B( Xo ,R) ^ — K> K > 0. Suppose u is any positive solution to 
the heat equation in B(xq,R) x [to — T, to]. Then 



\Vu\< 



ut 
u 



R 2 



c n i— |V«| 

+ — + c n K + c n VA sup 

T u 



c n |Vu| 

+ — sup 

R u 



in B(xq, R/2) x [t — T/2,t ]. Here c n depends only on the dimension n. 
Proof. 

By direct computation (see [H]), we have 



(A-ft)( 
In view of the estimate 



Vnl 



diudju 
OiOjU — 



diudju 

OoOjU — 

3 u 



> 



n 



u 



Au 



+ %Rij 



\Vu\< 



djudjii 



u 



u 



the above implies 



\Vu\< 



u 



> — | Au - 

nu 



|Vu| : 



+ 2Ri 



diudju 
u 
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Since An is also a solution to the heat equation, it follows that 

ts nw a |Vu|2 2/ |Vn| 2 \ 2 07 JVn| 2 
(A - d t )(-Au+ — ) > — ( An - - ) -2K- 

Let us write 



it nu\ u l u 



Then q satisfies 



Define 



Then H satisfies 



|Vn| 2 |Vn| 2 
-Au + = u t . 

u u 



(A-dt)q>— q 2 -2K^^. 

nu u 



H = q/u. 



(2.16) (A- d t )H >-H 2 -2K^- -2VHVlnu. 

n u 

Now we can use the Li-Yau idea of cut-off functions to derive the desired bound. The 
only place that may cause difficulty is that H may change sign. However it turns out 
that it does not hurt. Here is the detail. Let ip = 4>{ x ,t) De a smooth cut-off function 
supported in Qr^t = B(xq,R) x [to — T,to], satisfying the following properties 

(1) . V = il>(d(x,x ),t) = ip{r,t); ip{x,t) = 1 in Q r/2 ,t/4, < ip < 1. 

(2) . ip is decreasing as a radial function in the spatial variables. 

(3) . M < \m < g> when o < a < 1. 

(4) . |$ < 

Then, from (2.16) and a straight forward calculation, one has 
(2.17) 

>^.V(^) + 2^ U|2 

ip u 



A(ipH) - (ipH) t - 2-f ■ V(ipH) + 2ipK l —±- + 2V(^#)Vlnn 



2 Wib\ 2 
> -ipH 2 + (Aip)H - T—^H - ip t H + 2HVipV In u 
n ip 

= -4>H 2 - 2^^-H + (d 2 ip + (n - 1)— + d r M- log - V> t iT + 2tfVY>Vlnn. 

n ip r 

Suppose that at (y, s), the function ipH reaches a maximum. If the value is non-positive, 
there is nothing to prove. So we assume the maximum value is positive. Then (2.17) 
shows 



2ipK ] —^-+2 l —j^H > -ipH 2 +(d 2 ip+(n-l)^-+d r ipd r log x /7j)H-ip t H+2HVipV Inn. 



n z ip n r 

In the above, the only term we need extra care of is 

d r ipd r log y/gH. 
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Note that -C/R < d r ip/^ a < 0, d r log y/g < v 'K and H(y, s) > 0. Therefore 

12 IT7./.I2 _. v , 



|Vu| 2 |W| 2 ,- W 

2^i— J- + 2 l —^H + 2 v^i^-^V In i 



> — 4>H 2 + (d?iP + (n - l)^)H - CVk^H/R - ijj t H. 
nu r 



This shows that 



2 ,/|Vn| 2 u t \ 2 c n c n 2 ^|Vn| 2 ^c n , 2 

Hence 

|Vu| 2 c n c n T ^ /t-|Vu| c n |Vu| 

J — 2 1 " < 7§ + 7^ + ^ + L + ^ L. 

u^ u R z T u R u 

in the half parabolic cube. □ 

3. GRADIENT ESTIMATES ON THE LOG TEMPERATURE UNDER BACKWARD AND 

FORWARD RlCCI FLOW 

In this section we will prove certain localized or global gradient bound on the heat 
equation under backward and forward Ricci flow, i.e. equations (3.1) and (3.2) below. This 
estimate is a generalization of the results in [H] and [SZ], where the heat equation under 
a fixed metric is studied. Similar estimates for the conjugate heat equation (i.e. when A 
is replaced by A - R in (3.1) or (3.2)) were proven in [Ni3], [CKNT] and [CCGGIIKLLN] 
Chapter 8. This estimate then also relies on the derivative of the scalar curvature R. 

The current estimate under the forward Ricci flow ((3.2)) will be useful for Section 
5, where we will prove a global Gaussian upper estimate for Perelman's equation under 
nonnegative Ricci curvature assumption. 

Recall that the heat equation under backward and forward Ricci flow are given by 



(3.1) 

and 

(3.2) 

For (3.1) we have the following: 



Au - d t u = 0, 

An — dtu = 0, 
^tdij = —2Rij- 



Theorem 3.1. Let M be a compact Riemannian manifold equipped with a family of Rie- 
mannian metric evolving under the backward Ricci flow in (3.1). 
(a) (local estimate) . Suppose u is any positive solution to (3.1) in 

Qr,t = {(x, t) | x € M, d(x, x ,t) <R,t€[to- T, t ]} 

such that the Ricci > —k throughout. Suppose also u < M in Qr,t- Then there exists a 
dimensional constant c such that 

\Vu(x,t)\ . 1 1 f- , M v 

/ V < c(— + — T77T + V k) (l + log — A 

u(x,t) - { R T 1 / 2 n & u(x,t) J 

in Qr/2,t/2- 
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(b). (global estimate) Suppose u is any positive solution to (3.1) in M x [0, T]. Under 
the assumption that Ricci > 0, it holds 



\Vu(x,t)\ 1 , M 



U(x,t) ~ tl/2V 6 u(x,t) 

for M = sup Mx [ 0T ] u and (x,t) G M x [0, T]. 

Remark. As pointed out in [SZ], the local and global estimate can not replace each other. 
Also note that there is no other curvature assumption in part (b), nor any constants. 

Proof of Theorem 3.1 (a). 

We will use the idea in [SZ] with certain modifications to handle the changing nature of 
the metric. Suppose u is a solution to the heat equation in the statement of the theorem in 
the parabolic cube Qr,t- It is clear that the gradient estimate in Theorem 3.1 is invariant 
under the scaling u — ► u/M. Therefore, we can and do assume that < u < 1. 

Write 

/ = logu, w = |Vlog(l - f)\ 2 = . 

Since u is a solution to the heat equation, simple calculation shows that 

A/ + |V/| 2 -/ t = 0. 
We will derive an equation for w. First notice that 

_ 2V/(V/) f 2|V/| 2 / t 2Ric(Vf,Vf) 

Wt (I-/) 2 (I-/) 3 (I-/) 2 

_ 2V/V(A/ + | Vjf) 2| V/| 2 (A/ + | V/| 2 ) 2Ric(Vf, V/) 

(I-/) 2 (I-/) 3 (I-/) 2 

In local orthonormal system, this can be written as 

( j t_ (i-/) 2 + (i-/) 3 (i-/) 2 ' 

Here and below, we have adopted the convention fj = |V/| 2 and fa = Af. 
Next 



(3.4) Vw = (- 



+ 2- 



(l-/) 2 ^' (l-/) 2 (l-/)3- 
It follows that 



(3-5) = + t?^& + 



^f'ij , ^fifjjj ^fifijfj 
(I-/) 2 (I-/) 2 (I-/) 3 



^fifijfj Q /f/jj fl /? /j 



+ ,: ,v a + 2 ,; t '" +6 



(I-/) 3 (I-/) 3 (I-/) 4 ' 
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By (3.5) and (3.3), 

Aw — wt 

"(I-/) 2 (I"/) 2 

-I- fi LYZJ _j_ o fifijfj , cy_hjjj_ 
+ D (l-/)4 +5 (l-/)3 +/ (l-/)3 

_ ^ fifijfj _ n_f±fjj_ _ o l^/l 4 , ^Rjjfifj 
(I-/) 2 (I"/) 3 (I"/) 3 (I"/) 2 ' 

The 5th and 7th terms on the righthand side of this identity cancel each other. Also, by 
Bochner's identity 

fifijj ~ fjfiij = fjifjii ~ fiij) = Rijfifj- 

So the second term doubles with the last term. Therefore 
Aw — wt 

(3.6) 2/g |V/| 4 /i/p-Zi _ 4 iM. _ 2 |V/|4 + 

-(l_/)2+ (l-/)4 + (l_/)3 (1 _ /} 2 (l-/)3 + (l_/)2- 

Notice from (3.4) that 

2 fifijfj fifj 

Hence 

(3.7) = " 2V / Vw + 4 (|^3 ' 

(3.8) = -4^3 + PV/Vu, - 4 J^]^. 
Adding (3.6) with (3.7) and (3.8), we deduce 



Aw — wt 

2 f% , |V/| 4 , .fifijfj 



" (l- / )2 + (l-/)4 + (l_/)3 

+ -^V/W, - 2X7 fVw + 2 + l^M. 

^fij , n l^/l 4 , , fifijfj > n 
(I-/) 2 (I-/) 4 (I"/) 3 " ' 

2f iVfl 4 
A«i-w f > — ^V/Vw; + 2 . J ' - 4fcu;. 
1-/ (I-/) 3 
Since / < 0, it follows that 

2f IVfl 4 
Aw-w t > -^VfVw + 2(1 - f) ^_ f)4 ~ Mw, 



Since 
we have 
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i.e. 

2 f 

(3.9) Aw-w t > -^-j V/Vto + 2(1 - f)w 2 - Akw. 

From here, we will use a cut-off function to derive the desired bounds. Let tp = ip(x,t) 
be a smooth cut-off function supported in Qr,t, satisfying the following properties 

(1) . ^ = ij;(d(x, x ,t),t) = ijj(r, t); i/j(x, t) = 1 in Q r / 2 ,t/a, < ^ < 1. 

(2) . V is decreasing as a radial function in the spatial variables. 

(3) . M < ^ , J|36l < ^when < a < 1. 

< ?• 

Then, from (3.9) and a straight forward calculation, one has 
A(^jw) + b ■ V(ipw) - 2^- ■ V(V>u>) - (ipw)t 

> 2^(1 - f)w 2 + (b • Vip)w - 2^^-w + (Aip)w - ip t w - Aktpw, 
where we have written 

Comparing with the heat equation under a fixed metric, the last term —iptw is more 
complicated. It is given by 

dip , dip dd(x,x ,t) 

- Aw = - [ m + ^ at ] w - 

By our assumption that Ricci > — k and that — c/R < ^ < 0, we have 

—ip t w > ~~^ w ~ ckwtf) 1 / 2 . 

neie dt — dt . 
Therefore 

A(ipw) + 6 • V(V>u>) - 2—- • V(^w/) - (ipw) t 

(3.10) 2 

> 2^(1 - f)w 2 + (6 • Vip)w - 2^^-w + (AV>)w - - ckwip 1 / 2 . 

ip ot 

Suppose the maximum of ipw is reached at (x±,ti). By [LY], we can assume, without loss 
of generality that x\ is not in the cut-locus of M. Then at this point, one has, A(ipw) < 0, 
{ipw)t > and V(ipw) = 0. Therefore 

(3.11) 2%b(l-f)w 2 (x 1 ,t 1 ) < -[ (b-Vip)w-2^^-w+(Ajj;)w-^w](x 1 ,t 1 ) + ckw^ 2 . 
We need to find an upper bound for each term of the righthand side of (3.11). 
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|(6 . VVH < ^|V/|HW>| < 2™ 3 / 2 |/l ivvi 

/IW>| 



2[^(1 - f)w 2 f' A 



m ~ /)] 3/4 



<^(i_ /)u ,2 + (/IW-I) 4 



>(i-/)] 3 ' 

This implies 

f 4 

(3.12) | (b • WH < (1 - /)V^ 2 + c- 



'i? 4 (l-/)3- 

For the second term on the righthand side of (3.11), we proceed as follows 

IW| 2 /1/2 |Wf 

3.13 v , ,o 

/ 1 / 2^ JV^ 2 , 2 1 2 1 

Furthermore, by the properties of ip and the assumption of on the Ricci curvature, one 
has 

-(A^)w = -(d 2 ^ + (n - 1)^ + log ^5)w 



< (|a 2 V| + 2(n - + ^Vkw^ 

^ 1 , 2 . /rl 5 ^li2 , r 12 , ck 



Therefore 

(3.14) _(A^<^ 2 + c-^. 



Now we estimate \-^\ w. 



/1/2 1^ 
— w = ip ' tu- 



rn 



dt 1 r V 1/2 

1 iMi 

*s(* 1/ V'+4$)' 

This shows 

(3-15) |^-|u> < +cj^. 

Substituting (3.12)-(3.15) to the righthand side of (3.11), we deduce, 

2(1 - f)W < (1 - f)^w 2 + c R4{ (_ f)3 + + ^4 + ^2 + ^2 + ^V 7 ^- 
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Recall that / < 0, therefore the above implies 

/ 4 1 c c 

V>W 2 (xi,tl) < C 1 _ ^ + -^ 2 (xi,tl) + + ^ + Cfc 2 . 

f 4 

Since ntjji < 1 5 the above shows, for all in Qr,t, 

ip 2 (x,t)w 2 (x,t) < ip 2 (xi, ti)w 2 (x\, t\) 
< ip(x 1 ,t 1 )w 2 (x 1 ,t 1 ) 

C C ,o 

Notice that tp(x,t) = 1 in Qr/ 2 ,t/4 and w = |V/| 2 /(1 - f) 2 . We finally have 



We have completed the proof of Theorem 3.1 (a) since / = \og(u/M) with M scaled to 1. 



Proof of Theorem 3.1 (b). 

The proof is almost identical to that of Theorem 1.1 in [H] except for an additional 
curvature term. By direct computation, we have 



u u 



KJ I KJ j Li 



u 



In the above, comparing with the fixed curvature case (2.10), there is no more term 
containing the Ricci curvature. By (2.9), it holds 



A(wlog — ) - d t (u\og — ) = — k 

u u u 



Since 



(A-ft)(J^!!)> 



iVul' 



u 



the maximum principle implies that 

\Vu\ 2 



u 



1 M 

2 < T^log . 

t U 



u 



□ 



The remainder of the section deals with (3.2). For (3.2), we no longer have the nice can- 
celation effect that associated with (3.1). So we only obtain the following global gradient 
estimate under curvature assumptions. 

Theorem 3.2. Let M be a complete Riemannian manifold equipped with a family of 
Riemannian metric evolving under the forward Ricci flow in (3.2) with t 6 [0, T]. Suppose 
u is any positive solution to (3.2) in M x [0, T]. Then, it holds 

\Vu(x,t)\ < fl I M 
u(x,t) - V t Y ° g u(x,t) 

for M = sup Mx [ 0T ] u and (x,t) £ M x [0, T]. 
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Moreover, the following interpolation inequality holds for any 5 > 0, x, y G M and 
< t < T: 

u(y,t) < c 1 n(x,t) 1 /( 1 + 5 )M 5 /( 1+<5 )e C2d ^ 2 / t . 
Here c\,C2 are positive constants depending only on 5. 

Proof of Theorem 3.2. 

This again is almost the same as that of Theorem 1.1 in [H] . By direct calculation 

A(itlog — )-d t (ulog— ) = -■ ! 



(A- 8 «)( |V " 12 - 2 



u u u 

2 

> 0. 



u u 



diudju 

OiOjU — 

J u 



The first inequality follows immediately from the maximum principle since 

\Vu\ 2 M 
t- ■ it log — 

u u 

is a sub-solution of the heat equation. 
To prove the second inequality, we set 

l(x,t) = \og{M/u(x,t)). 

Then the first inequality implies 

\Vy/l(x,t)\ < 1/Vt. 
Fixing two points x and y, we can integrate along a geodesic to reach 

Vlog(M/n(x,t)) < y/\og(M/u(y,t)) + ^E^l. 

The result follows by squaring both sides. □ 

4. POINTWISE AND GRADIENT ESTIMATE FOR THE FUNDAMENTAL SOLUTION TO AN 

EQUATION OF PERELMAN'S 

In the paper [P] Perelman introduced an equation which after time reversal becomes 

J Ait — Ru — dtu = 0, 



(4.0) 



dtdij — ^Rij- 



Here as before A is the Laplace-Beltrami operator with respect to the metric g^j evolving 
by the backward Ricci flow. R is the scalar curvature. This equation and the associated 
monotonicity formula have proven to be of fundamental importance. Using the maximum 
principle and reduced distance, Perelman proved a lower bound for the fundamental so- 
lution to (4.0). An outstanding feature of the estimate is that it does need any explicit 
curvature assumption. The information on curvature is encoded in the reduced distance. 
From the analysis point of view, it would be desirable to establish an upper bound for the 
fundamental solution too. Here we first prove an upper bound under no explicit curvature 
assumptions. The bound is in terms of more traditional geometric quantities, i.e. the best 
constant in Sobolev imbedding or Yamabe constant, which are controlled by the lower 
bound of the Ricci curvature and injectivity radius. Under more restrictive curvature as- 
sumptions, we are able to prove a Gaussian like upper bound. Let us mention the method 
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by maximum principle alone does not seem to yield the upper bound. Our method is 
based on the one by J. Nash. For related results on local lower and upper bounds for 
fundamental solutions of (3.2) and for a global lower bound for the conjugate of (3.2) in 
the spirit of Perelman, please see the interesting papers [G], [Nil] and [Ni2]. 

In order to state our theorem, we need to recall two concepts. One is the Yamabe 
constant and the other is the best constant in the Sobolev imbedding. 

Given a Riemannian metric g{t) the Yamabe constant is 

Y(t) = mj -. — ^-7—. 

(/0V(n-2) d/i ( M ))("- 2 )/" 

The other is a Sobolev imbedding theorem due to E. Hebey Heb] which is a refined form 
(on the controlling constants) of the result by T. Aubin ACDH : 

Theorem S. Let M be a complete (compact or noncompact) Riemannian n-manifold. 
Suppose the Ricci curvature is bounded below by k and the injectivity radius is bounded 
below by i > 0. For any e > 0, there exists B(g) = B(e,n, k,i) such that for any eft £ 

r \ (n-2)/n 

\ \u\ 2n l^d^g)) <(K(n) 2 + e) \Vu\ 2 dfi(g) + B(g) u 2 dfi(g). 

Jm / Jm Jm 

Here K{n) is the best constant in the Sobolev imbedding in R n . 

We also need to mention the result by Hebey and Vaugon [HV] where Theorem S is 
proven with e = 0. However, then the constant B may depend on the derivative of the 
curvature tensor which is harder to control. For our purpose, it suffices to fix the e as any 
positive constant, say 1. 

The following theorem is the main result of the section. It contains three statements. 
The first one is an upper bound controlled by the Yamabe constants, the second is an 
upper bound controlled by the constant B(g) in the Sobolev imbedding Theorem S. They 
may seem technical at the first glance. However, the third statement of the theorem 
provides a clarification. It shows that these upper-bounds are the proper extension of 
on-diagonal upper bound for the heat kernel in the fixed metric case. Recall that for 
a compact Riemannian manifold M without boundary, the heat kernel G satisfies the 
following on-diagonal upper bound: 

G(x,t;y,s) < c x max{ Tj—K^ A} 

for some constant c%, ci > and for all t > s and x, y £ M. 

Here are some additional notations for the theorem. We will write 

R = -mm{R(x,t),0} 

where R(x,t) is the scalar curvature under the metric g(t). When the scalar curvature 
changes sign, the theorem will also involve the expression 

1 



(maxi?-(-,t))- 1 + (t- s)' 
This quantity is regarded as when R(-,t) > 0. 
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Remark 4.1. In statements (b) and (c) of Theorem 4.1 below, the controlling constants 
depend only on the dimension, the lower bound of Ricci curvature and the lower bound 
of injectivity radii. By the result of Cheeger [Che], if one assumes that the sectional 
curvatures are bounded between two constants and the volume of geodesic balls of radius 
1 is bounded below by a positive constant, then the injectivity radii are bounded from 
below by a positive constant. Therefore, the controlling constants in (b) and (c) depend 
only on the bound of sectional curvature, the lower bound of volume of balls of radius 
1 and dimension. The same can be said for Theorem 5.1 below. The upshot is that the 
length of time and the incompatibility of metric at different time do not destroy the bound. 

Theorem 4.1. Suppose equation (4-0) has a smooth solution in the time interval [s,t] 
and let G be the fundamental solution of (4-0). Then the following statements hold. 

(a) . Suppose the Yamabe invariant Y(g(r)) > for r £ [s,t], then 

G(x,t;y,s) 

— / \ n/4 ' 

/ j(t+s)/2 e2Cn a(T)/n Y ( T )dT /* +s)/2 e -^ n a(r)/n [l + ( t _ T ) max R~ (-, t)] ~ 4/ " Y{r)dr J 

Here a(r) = £ {maxR - { .i))^+{t-i) dL 

(b) . Let B(g(r)) be the best constant in the Sobolev imbedding Theorem S. Then 

G(x,t;y,s)< 



^j(t+s)/2 e 2H(r)c n /n dT / ( i t+s)/2 [l + (t- r ) maXj R-(.,t)]-4/n e -2H(r) Cn /n dr ^ 



n/4 



with 



(c). In the special case that R(-,t) > and Ric(g(r)) > k and the injectivity radius is 
bounded below by i > 0, for all r £ [s,t], then 

G(x, t; y, s) < C(n, B) max{ - _^ n/2 , !}■ 

Here B only depends on n, k and i. Moreover 

\VyG(x,t;y,s)f 1 ^j^j 

G(x,t;y,s)* ~ 1 ' ] (t - s) S G(x,t;y,s) ' 

Remark 4-2. Recently, in a paper [CL], Chang and Lu, proved a derivative estimate for 
the Yamabe constant under the Ricci flow. It can be coupled with this theorem to obtain 
better upper bound on G. 

Professor Lei Ni also informs us that he also knows a result on upper bound in the case 
of certain Sobolev inequality. 

Proof of part (a). Without loss of generality, we take s = here and later. 
Let G be the fundamental solution to (4.0). By the reproducing property 

G(x,t;y,0) = J G(x,t;z,t/2)G(z,t/2;y,0)dn(z,t/2), 
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there holds 
(4.1) G(x,t;y,0)< 



1/2 



1/2 



G 2 (x,t;z,t/2)dfi(z,t/2) J G 2 (z,t/2;y,0)dfi(z,t/2) 

Therefore an upper bound follows from pointwise estimate on the two quantities 

(4.2) p(t) = J G 2 (x,t;y,s)d^x,t), 

(4.3) q(s) = J G 2 (x,t;y,s)dfi(y,s). 

Let us estimate p{t) in (4.2) first. 
It is clear that 

j t p(t) = 2 j G[AG — RG]dfx(x,t) + J G 2 Rdn{x,t). 

Here and later we omit the arguments on G and differential operators when no confusions 
appear. Therefore 

j t p(t)<-j[\VG\ 2 + RG 2 W(x,t) 

(^•4) p n — 2 3n — 2 f 

= -j[\^G\ 2 + ^—^RG 2 ]d^x,t) - ^— ^ J RG 2 d^x,t). 

Let Y(t) be the Yamabe constant with respect to g(t), i.e. 

. jm\ 2 + ^R4> 2 W{x,t) 

Y(t) = mf -, — 7^i—. 

(J<f>^-Vdv(x,t)) {n - 2)/n 

By Holder's inequality 

J G 2 d^x,t) < [ j G 2n ^ n - 2 Uf,(x,t)] {n - 2)/{n+2) [ j G^(x,t)] 4/(n+2) , 

we arrive at the 'conformal' Nash inequality 
(4.5) 

J G 2 dfi(x,t) 

<c n Y{t)-^) [|[|V0| 2 + ^-^^ 2 ]^(x,t)] n/( " +2) [ / GdKx,t)] 4/{n+2) . 
It is easy to check that 

J G(x,t;y,s)dn(x,t) = 1. 

Therefore (4.5) becomes 

J G 2 d^(x,t)c n < CnYit)-"^ [|[|V0| 2 + _^L/^]d M (x,t)] B/(B+2) . 

Substituting this to (4.4), we deduce 

p'(t) < - Cn p(t) (n+2)/n Y(t) - [ RG 2 d^(x,t). 

4(n -I) J 
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It well-known (see [CK] e.g.) that the scalar curvature R satisfies the inequality 
This implies 



^ + A J R + c n i? 2 <0. 
at 



R{y,T) ~ (maxi?-(-,t))-i + c n (t - r)' t <C t. 
Here and later, if R(-,t) > 0, then the above fraction is regarded as zero. Hence, for 
r G (s,t), 

*Y < -«,(r)<«>/«r(r) + 1 — jF ^L ZT - 1 — f(T) . 

Let 

r i 

a(r) = / — — — — — j — dl. 



(maxR-(-J))- 1 + (t-l) 
Then the above ordinary differential inequality becomes 

(e- c "°( T )p(T))' < c n (e- c " a(T V(r)) (n+2)/n e 2c " a(T)/n y(r). 
Integrating from s to t, we deduce 



[ J s *e 2c « a «/"y(r)dr]" /2 
This immediately shows that 



(4.6) y G 2 (x, t; y, s)^(x, t) = p(i) < 



CnC c n a(t) 



[jt e 2c n a(T)/n Y ( T )d T 



n/2 ' 



Next we estimate q(s) in (4.3). Due to the asymmetry of the equation, the computation 
is different. Notice that the second entries of G satisfies the backward heat equation, i.e. 

A y G(x, t; y, s) + d s G(x, t; y, s) = 0. 

This gives 

q'(s) = -2 J GAGdv(y,s) + J RG 2 d f i(y,s). 

Hence 

(4.7) q'{s)> j[\VG\ 2 + RG 2 W(y,s). 

By the same argument as before we arrive at the Nash inequality 
(4.8) 

<c n Y{s)-^ 2 ) [ J [m 2 + ^^R^ 2 W(y,s)] n/in+2) [ J Gd^y,s)] 4/{n+2 \ 
This time we have to compute the quantity 

I(s) = J G(x,t;y,s)dfj,(y,s). 
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It is clear that 

(4.9) I'(s) = J G(x, t- y, s)R(y, s)d^y, s). 
Recall that 

R(y,T) > — — — ; — — — ; -, t <t. 

yy > 1 - (max i?- (•,*))-! + Cn (t- r )' 

Combining this with (4.9) we deduce 

1 M " "(max^(-,t))- 1 +c n (t-r) /(T) ' 
Integrating from s to t and noting that I(t) = 1, we obtain 

(4.10) I(s)<l + Cn(t-s)maxR~(;t). 

Substituting (4.10) to (4.8), we deduce 
(4.11) 

[\VG\ 2 + RG 2 W(y,s) 

/n — 1 In — 2 f 

[\VG\ 2 + ^—^R&] M y, s) + _ - J RG 2 dfi(y, s) 

> [g( s )](«+ 2 )/"[l + Cri (t- s )ma X i?-(.,t)]- 4/n Y(s) Cnq{s) 



(max J R-(-,t))- 1 + (t- s)' 

Here, again we used the lower bound on the scalar curvature, given just below (4.9). 
This and (4.7) together imply that, for r € (s,t), 

«r) > * W r)](-«W.[l + ( t -r)m„ir(,f)]-^ y(r)- ^-(.^ + (t _ r) «M- 



Let again 



(max J R-(-,t))- 1 + (t-/) 
Then 

(n+2)/n 



g(r)e 



c„a(r) 



g(r)e c " a{r) e" 2c " a(T)/n [l + (t - r) maxir(-, t)] ~ 4/ ™ y(r). 



Integrating from s to i, we obtain 

q(s) = J G 2 (x,t;y,s)dfi(y, s ) 
(4.12) c ra e- c " Q ( s ) 



^ £ e -2c n a(r)/n [l + (* _ T ) max i?" (•,*)] ~ 4/ ™ F(r)(iT^ 

Now (4.6) an (4.12) respectively imply that 

/„ „c n a(t/2) 
G 2 (z,t/2;y,0)dfi(z,t/2) < 



n/2 ' 



Jl' 2 e 2c ^/ n Y(T)dT 
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G 2 (x,t;z,t/2)dfi(z,t/2) 



(4.14) Crae -c n aft/2) 



( /t/2 e- 2c ^)/n [i + (t- r) max i?" (■,*)] _4/n F (r)dr^ 



n/2' 



n/4" 



By (4.1), (4.13) and (4.14), we arrive at the following upper bound 

G(x,t;y,0) 

^ J * /2 e 2c " a (-)/™y(r)dT // /2 e -2c n a(r)/n [i + (t _ T ) maxi?-(-, *)] _4/n Y(r)^ 

This proves part (a). 
Proof of Part (b). 

We generally follow the previous arguments between (4.1) and (4.14) to derive an upper 
bound. The difference is that we will use the Sobolev inequality (Theorem S) instead of 
the Yamabe constant. 

As before, by Holder's inequality and Theorem S, we arrive at the Nash type inequality 

f G 2 dn{x,t) 
(4.15) J 

< [ / [c n \VG\ 2 dfi(x,t) + B(g(t)) / G 2 dfi(x,t)] n/{n+2) [ / Gdfi(x,t)} 4/{n+2) . 
J Jm J 

Here and later c n is a dimensional constant that may change from line and to line. Since, 
again, 

1, 



J Gdn(x,t) 



we have 

/- -i (n+2)/n 

\VG\ 2 dn(x,t)>c n / G 2 dn(x,t) - c n B(g{t)) G 2 dfj,(x,t). 

L JM J J 

Combing (4.16) with (4.4) under again the notation (4.2), we obtain 

At) < -c nP (t) {n+2)/n + c n B(g(t))p(t) - [ RG 2 dn(x,t). 

Jm 

Fixing s and t, for any r G (s,t), we still have the lower bound for the scalar curvature 
(just after (4.9)) 

R(-,t) > —- — — r^— : ; -, T <t. 

(maxi?-(-,t))- 1 + c n (t-r)' 

Therefore 

(4.17) p\r) < -c n p(r) (n+2)/n + Cnh(T) P (r), 

with 

h(r) = B(g(r)) + ' 



(max R~ (•, t))- 1 + (t-r)' 
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Let H(t) be the anti-derivative of h(r) such that H(s) = 0. Then 

(n+2)/n 



" c " I(T) p(r)J < -c n (e- c ^Mp(r)y 



e 2c n H{r)/n_ 



Integrating from s to t, we arrive at 

p c„ff(t) 

(4.18) p(t) < 



n/2 ' 

J* e 2H(r)c n /n dT ' 



Our next task is to bound 



q(s) = J G 2 (x,t;y,s)dfi(y,s). 
Clearly the counter-parts of (4.7), (4.10) and (4.15) still hold. i.e. 

q'(s)> j(\VG\ 2 +RG 2 )d^y,s). 

I{s) = j Gdfi(y, s) < 1 + c n (t — s) maxi?~(-, t). 

j G 2 d^y,s) 

< [ j [c n \VG\ 2 d^y,s)+B(g(s)) j G 2 d^x,t)} n/{n+2) [ j Gd^y, s)] 4/(n+2) 

Also 



R(y,s) > — , 



[maxi?~(-, i)] _1 + c n (t — s) ' 
These four inequalities imply that 

q'(s) > c n q(si n+2 ^ n [l + {t-s)^R'{;t)Y A l n -c n h{s)q{s). 

Here h(s) is given by the expression just below (4.17) with r replaced by s. Now, for fixed 
s and t and any r £ (s,i), the above differential inequality on q'(s) is still valid for q'{r) 
when s is replaced by r. Let H(t) be the antiderivative of h{r) with H(s) = 0. Then it 
is clear that 

/ V / \ ( n+2 )/ n 

(e c " H ^q(r) J > c n (e c " H ^q(r)\ [1 + (t - r) max iT(-, t )]~ 4 / n e - 2H ^ c "/ n . 

Integrating from s to t, we arrive at 

r p-CnH(s) 

*W < 7 W^- 

( /«[! + (t-r)maxi?-(-,t)]" 4 /n e -2//(r)c„/n dr j 

By (4.18), we have 

Cne c^/2) 



(4.20) p(*/2) = / G 2 (z,t/2;y,0)^(z,t/2) < — - 



n/2 ' 

e 2//(r)c„/n dr ' 
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Also, (4.19) shows 
(4.21) 

Cn.C 



q (t/2) = j G 2 (x,t;z,t/2)dfi(z,t/2) < 



Here 

ft/2 



Jj 2 [l + (t- T)maxK-(.,t)]- 4 /» e ~2%)/n ( | r j 



n/2' 



ff«/2)-_( [%Mi t ( „ ai ,;, t „_ r| ^ 

Multiplying (4.20) and (4.21), and using (4.1), we have proven the on-diagonal upper 
bound 

G(M;y,0) 2 < -75- ^ -j- 2 

jt/i e 2H(r)c n /n dT \ t jt^ + (t - T ) max R~ (•, t)}-^ n e~ 2H ^ l n d T J 

This gives part (b). 
Proof of part (c). 

In the special case that R(-,t) > and Rice(g{r)) > fc uniformly and the injectivity 
radius is uniformly bounded below by i, then 

Htr, = JjB W) ) + (m „ fl - (%t ; ) - 1 + (t _ i) H - ^(»,M)r 

and 

1 



(max J R-(-,t))- 1 + (t - t) 
Hence the above immediately shows 

G(x,t;y,0) <C(c n , B) max{^,l}. 

The gradient estimate follows from Hamilton's argument in Theorem 1.1 [H], which can 
be easily generalized to the present case. For z € M and r 6 [0, t/2], let 

r) = G(x, t; z, r). 

Then v is a solution to the backward heat equation At; + v T = 0. By direct computation 



<A + *)<^)-? 
> 0. 



ajajt; — 



Let ^4 be the maximum of v in the time interval [0, t/2]. By the above estimate 

A < C(c n , 5) max{^,l}. 
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Direct computation shows 



A A 
A(v log — ) + d T (v log — ) 
v v 

= (Av + d T v) log A + (A + d T )(v log—) 

\Vv\ 2 

v 

Let 4> = (t/2) — t, then it is clear that, for 

ft = flog — , 

there holds 

Aft, + <9 r /i > 0. 

By the maximum principle, applied backward in time, we have 



< C- log 



V 2 T V 

for r £ [0, t/4]. □ 

5. The case of nonnegative Ricci curvature 

In this section, we specialize to the case of nonnegative Ricci curvature. We establish 
certain Gaussian type upper bound for the fundamental solution of (4.0). We will begin 
with the traditional method of establishing a mean value inequality via Moser's iteration 
and a weighted estimate in the spirit of Davies [Da]. However, there is some difficulty 
in applying this method directly due to the lack of control of the time derivative of the 
distance function. The new idea to overcome this difficulty is to use the interpolation 
result of Theorem 3.2 and the bound in Theorem 4.1 (c). 

The following are some additional notations for this section. We will use B(x,r;t) to 
denote the geodesic ball centered at x with radius r under the metric g(t); \B(x,r;t)\ s to 
denote the volume of B(x, r; t) under the metric g(s). 

The main result of this section is Theorem 5.1 below. Note that the theorem is quali- 
tatively sharp in general since it matched the well-known Gaussian upper bound for the 
fixed metric case. Also there is no assumption on the comparability of metrics at different 
times. In this theorem, we assume the manifold is compact. This accounts for the extra 1 
on the Gaussian upper bound. Even in the case of fixed metric, the heat kernel converges 
to a positive constant for large time. The theorem still holds for certain noncompact 
manifolds under suitable assumptions. In this case the extra 1 in the upper bound should 
be replaced by 0. 

Remark 5.1. As mentioned in section 4 (Remark 4.1), the controlling constants in the 
theorem below can be made to depend only on the bound of sectional curvature, the lower 
bound of volume of balls of radius 1 and the dimension. 

In the case Ricci > —k with k > 0, then certain integral Gaussian bound similar to 
the one below (5.13) can still be proven by the same method. However, so far we are not 
able to derive a pointwise Gaussian upper bound without an exponentially growing term 
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e kt . This is due to a lack of an efficient mean value inequality for the second entries of the 
fundamental solution, which satisfies (3.2) after a time reversal. 

Theorem 5.1. Assume that equation (4-0) has a smooth solution in the time interval 
[0, T] and let G be the fundamental solution of (4-0). Suppose that Ricci > and that 
the injectivity radius is bounded from below by a positive constant i throughout. Then the 
following statement holds. 

For any s,t € (0, T) and x,y £ M, there exist a dimensional constant c n , a dimension 
less constant c and a constant A depending only on i such that 

Proof 

It is obvious that we only have to deal with the case that B(x, 2^t — s, s) is a proper 
sub-domain of M. Otherwise, y/t — s/2 > d(x,y,s) for any x,y £ M. So the exponential 
term is mute and the result is already proven by Theorem 4.1 (c). 

First we use Moser's iteration to prove a mean value inequality. The only new factor 
is a cancelation effect induced by the backward Ricci flow. So we will be brief in the 
presentation at this part of the proof. 

Let u be a positive solution to (4.0) in the region 

Qcrr{x, t) = {(y, s) \ z e M, t — (ar) 2 < s < t, d(y, x, s) < ar}. 
Here r > 0, 2 > a > 1. Given any p > 1, it is clear that 
(5.1) Au p - P Ru p - d t u p > 0. 

Let 4> : [0, oo) — ► [0, 1] be a smooth function such that \(j>'\ < 2/((<j— l)r), ft < 0, <f> > 0, 
4>(p) = 1 when < p < r, (p(p) = when p > ar. Let r\ : [0, oo) — ► [0, 1] be a smooth 
function such that |r/| < 2/((cr — l)r) 2 , > 0, r] > 0, i^(s) = 1 when t — r 2 < s < t, 
(j)(s) = when s < t — (ar) 2 . 

Writing w = u p and using wip 2 as a test function on (5.2), we deduce 

(5.3) J V(wip 2 )Vwdp(y, s)ds + p J Rw 2 ip 2 dp(y, s)ds < — J (d s w)wip 2 dp(y, s)ds. 
By direct calculation 

(5.4) J V(wtp 2 )Vwdp(y,s)ds = J \V(wip)\ 2 dp(y, s)ds - J \Vip\ 2 w 2 dp(y, s)ds. 
Next we estimate the righthand side of (5.3). Here we will use the backward Ricci flow. 
-J (d s w)wip 2 dp(y,s)ds = J w 2 vbd s ipdp(y, s)ds+^ J (wifj) 2 Rdp(y,s)ds-^ J (wifj) 2 dp(y,t). 
Observe that 

d s ?p = r](s)(t)'(d(y, x, s))d s d(y, x, s) + (/){d(y, x, s))rj'(s) < 4>(d(y, x, s))rj'(s). 
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This is so because <f>' < and d s d(y, x, s) > under the backward Ricci flow with nonneg- 

ative Ricci curvature. Hence 

(5.5) 

(d s w)wip 2 d/j,(y , s)ds 



< J w 2 ^<j)(d(y,x,s))r]'(s)dfi(y,s)ds + ^ J ' (wip) 2 Rdfj,(y, s)ds - ^ J " (wi/j) 2 d/j,(y , t) . 
Combing (5.3) to (5.5), we obtain, in view of p > 1 and R > 0, 

(5.6) / |V(wY>)| 2 <f/x(y,s)cfe + - I (wi/j) 2 dn(y,t) < _° 2 / w 2 dfi(y, s)ds. 
J z J [a i) r jQ ar(x t) 

By Holder's inequality 
(5.7) 

// /• \ (n-2)/n / \ 2/n 

(^)2(l+(2/n)^ (y) s) < / y ty„,)V(»-2)) dM(y) s) j ^ y (^) 2 d/x(y, S )J . 

Let us assume that B(x, ar, s) is a proper sub-domain of M. In this case, for manifolds 
with nonnegative Ricci curvature, it is well-known that the following Sobolev imbedding 
holds (see [Sa] e.g.) 

/\ (n— 2)/n 2 2 r 

^ w fn/(n-2) My ^ s) < C " CTr 2 . / [| V0M| 2 + r" 2 (^) 2 ]^(y, S). 

/ \B(x,ar,s)\ s ' J 

For se[t- (crr) 2 ,t], by the assumption that the Ricci curvature is nonnegative, it holds 

B{x,ar,s) D B(x,ar,t); \B(x, ar, s)\ s > \B(x,ar,t)\ t _^ ar ^2. 

Therefore we have 
(5.8) 

/\ (n— 2)/n 2 2 r 

/ |S(x,(rr,t)| t i K)2 ^ 

For s G [i — (crr) 2 ,i]. Substituting (5.7) and (5.8) to (5.6), we arrive at the estimate 

[ w 29 dn(y,s)ds < c n — ^ (- / w 2 d[i(y,s)ds 

jQr{x,t) \B(x,ar,t)\ 2/ ™ (ar)2 V( CT ~ i)^ 2 ^Q CT ,(x,t) 

with 6* = 1 + (2/n). Now we apply the above inequality with the parameters o"o = 1, <?{ = 
1 — S* =1 2~ J, ~ 1 and p = 9 l . This shows a L 2 mean value inequality 

(5.9) sup u 2 < Cn / u 2 d^(y,s)ds. 

Q r/2 (x,t) r 2 \B{x,r,t)\ t _ r 2 J Qr{x , t ) 

From here, by a generic trick of Li and Schoen [LS], applicable here since it uses only the 
doubling property of the metric balls, we arrive at the L 1 mean value inequality 

(5.10) sup u < ° n / udn(z,T)dT. 

Q r/2 (x,t) r 2 \B{x,r,t)\ t _ r 2 J Qr{X:t) 
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Fixing y G M and s < t, we apply (5.10) on u = G(-, •; y, s) with r = y/t — s/2. Note 
that J M u(z,T)d/j,(z,T) = 1. The doubling property of the geodesic balls show that 

G(x,t;y,s) < - — = — — 

\B{X, y/t - S,t)\ s 

when \B(x, y/t — s, s)\ is a proper subdomain of M. 

Without loss of generality, we take s = 0. We begin by using a modified version of the 
exponential weight method due to Davies [Da]. Pick a point xq € M, a number A < and 
a function / G L 2 (M, <?(0)). Consider the functions F and u defined by 

F(x,t) = e Xd ^' X0 ^u(x,t) = e ^x,x ,t) J G(x,t;y,0)e- Xd ^ xo ^f( y y)dii( y y,0). 
It is clear that u is a solution of (4.0). By direct computation, we have 

= 2Ay e 2Xd ^ x ' Xo ^d t d(x,x ,t)u 2 (x,t)dn(x,t) + j e 2Xd ^ x °^u 2 (x,t)R(x,t)dn(x,t) 

+ 2 y e 2Ad ( x ' X0 '*)[Au-i2(x,t)«(x,*)]u(x,*)d/i(»,*)- 
By the assumption that Ricci > and A < 0, the above shows 

d t j F 2 (x,t)dfi(x,t) < 2 j e 2Xd ^ X0 ^uAudfi(x,t). 
Using integration by parts, we turn the above inequality into 
d t J F 2 (x,t)dfi(x,t) 

< -4A J e 2Xd ^ x ' X0 ^uVd(x,x ,t)Vudn(x,t)-2 j e 2Xd( - x ' X(ht) \Vu\ 2 dfi(x,t). 
Observe also 

J\VF{x,t)\ 2 d^{x,t) = J \V(e Xd ^ x ' X0 ^u(x,t))\ 2 dfi(x,t) 

= J e 2Xd ^ x ' X0 ^\Vu\ 2 dn(x,t) + 2X J e 2Xd( - x ' X0 ^uVd(x,x ,t)Vudn(x,t) 

Combining the last two expressions, we deduce 

d t j F 2 (x,t)dfi{x,t) < -2 J \VF{x,t)\ 2 dfi{x,t) + X 2 J e 2Xd ^ x ' X0 ^\Vd\ 2 u 2 dn{x,t). 
By the definition of F and u, this shows 

d t J F 2 (x,t)dfi(x,t) < X 2 j F(x,t) 2 dfi(x,t). 
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Upon integration, we derive the following I? estimate 

(5.11) / F 2 (x,t)dfi(x,t) < e xh \ F 2 (x,0)dn(x,0) = e xh \ f (x) 2 dfi(x , 0) . 



Recall that u is a solution to (4.0). Therefore, by the mean value inequality (5.9), the 
following holds 

u(x,t) 2 < p — / / _ u 2 (z,r)dfi(z,T)dT. 



t\B( X , t)\ t/2 Jt/2 JB(x,y/t/2,T) 

i.e. By the definition of F and u, it follows that 



u(x,t) 2 < C " r f e- 2 ^^F 2 {z,r)dn{z,r)dr. 

t\B(x, Wt 2, t)\ t /o Jt/2 JB(x^/tj2,T) 



t\B(x,y/t/2,t)\ t /2 Jt/2 J B(x,yftj2,r) 

In particular, this holds for x = xq. In this case, for z € B(xq, y/t/2, r), there holds 
d(z,xo,r) < \Jt/2. Therefore, by the assumption that A < 0, 



e -2\y/tj2 nt 

u(x ,t) 2 < = / / _ F 2 (z,T)dfi(z,T)dT. 



t\B(x ,y/t/2,t)\ t /2 Jt/2 JB{x ,y/t/2,T) 

This combined with (5.11) shows that 

\ 2 t-\V2l r 

u( Xo ,t) 2 < n / m 2 Mv,o) 

\B(x , \/t 2,t ) t/o J 



\B(x ,y/t/2,t)\ t/2 
i.e. 

(5.12) G(x ,t-,z,0)e- xd ^^f(z)d f i(z,0)) < n / f(y) 2 d»(y,0). 

\J J \B(x , y/t/2,t)\ t / 2 J 

Now, we fix yo such that d(yo,xo,0) 2 > 4a 2 t with a > 1 to be chosen later. Then it is 
clear that, by A < and the triangle inequality, 

-Xd(z,x o ,0) > -aXd(x o ,y o ,0) 

when d(z,yo,0) < \ft. In this case, (5.12) implies 
(5.13) 

/ G(x ,t;z,0)f(z)d^z,0)) < JL — /=— f(y) 2 d^y,0). 

\JB(y ,Vi,0) J \B(x ,^/t/2,t)\ t/2 J 

Now we take 

__ d(x o ,y o ,0) 
X ~ bt • 

Take b > and a > sufficiently large. Then (5.13) shows, for some c > 0, 



/ 



G 2 (x ,t;z,0)d(i(z,0) < 



„ „-cd(xo,yo,0) 2 /t 



HvcVtfi) \B(x ,y/t/2,t)\ t /2 
Hence, there exists zq G B(yo, \/t, 0) such that 



G 2 (x ,t;z ,0) < 



-cd(x ,y ,0y/t 



B(x ,y/i/2,t)\ t/2 \\B(x ,y/i,0)\ 
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By the doubling property of the geodesic balls, it implies 

-cd(x ,y ,0) 2 /t 

\B(x ,Vt,t)\o\\B{x Q ,Vt,0)\o 

Finally, let us remind ourself that G(xq, t; - , ■) is a solution to the conjugate equation of 
(4.0). i.e. 

A z G(x, t; z; r) + d T G(x, t; z, r) = 0. 

Therefore Theorem 3.2 can be applied to it after a reversal in time. Consequently, for 
5>0,C>0, 

(5.15) G(x ,t; y , 0) < CG 1 '^ (x , t, z , 0)M S /^ , 

where M = sup Mx ^ G(xo, t, ■, ■). By Theorem 4.1, part (c), there exists a constant 
A > 0, depending only on the lower bound of the injectivity radius such that 

M < ^maxj^TT, 1}. 

This, (5.14) and (5.15) show, with 6=1, that 

C A e -cd(x ,y ,0) 2 /t 

G(x Q , t; y , 0) 2 < max{— ^, 1}- 



^" /2 ' ^\B{x ,Vi,t)\o\\B(xo,Vt, 0)| 
By the assumption that the Ricci curvature is nonnegative, we have 

\B(x ,Vi,t)\o < \B(x ,Vi,0)\o. 

Therefore 

C A e -cd(x ,y ,0) 2 /t 

GW yo ,0)<max { ^,l} — ^ - . 

Consequently 

G (,„, iiw ,0, < ^(l + ^ + ^_l^)e-<»-».^. 
Since xo and yo are arbitrary, the proof is done. 

□ 
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